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ABSTRACT 

Kouveliotou et al. (1993) recently confirmed that gamma-ray bursts are 
bimodal in duration. In this paper we compute the statistical properties 
of the short (< 2 s) and long (> 2 s) bursts using a method of analysis 
that makes no assumption regarding the location of the bursts, whether in 
the Galaxy or at a cosmological distance. We find the 64 ms channel on 
BATSE to be more sensitive to short bursts and the 1024 ms channel is more 
sensitive to long bursts. We show that all the currently available data are 
consistent with the simple hypothesis that both short and long bursts have 
the same spatial distribution and that within each population the sources 
are standard candles. The rate of short bursts is ~ 0.4 of the rate of long 
bursts. Although the durations of short and long gamma-ray bursts span 
several orders of magnitude and the total energy of a typical short burst is 
smaller than that of a typical long burst by a factor of ~ 20, surprisingly the 
peak luminosities of the two kinds of bursts are equal to within a factor of 
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1 Introduction 



The Burst and Transient Source Experiment (BATSE) on board the Comp- 
ton Gamma-Ray Observatory has demonstrated that the distribution of 
gamma-ray bursts is isotropic over the sky and bound in the radial direc- 
tion (Meegan et al. 1992). This strongly indicates that the sources of the 
bursts are located either at cosmological distances (Paczyhski 1991, Dermer 
1992, Mao & Paczyhski 1992, Piran 1992) or in an extended Galactic halo 
(Li & Dermer 1992, Duncan, Li, & Thompson 1992). Many models have 
been proposed for the cosmological scenario (e.g., Eichler, et al., 1989, Pi- 
ran, Narayan, & Shemi, 1992, Narayan, Paczyhski, & Piran 1992, Rees & 
Meszaros 1992, Usov 1992, Woosley 1993) as well as for the extended halo 
scenario (Fabian & Podsiadlowski 1993). For a recent review on gamma-ray 
bursts, see Paczyhski (1992). 

Recently Kouveliotou et al. (1993) showed that the distribution of dura- 
tions of gamma-ray bursts is bimodal. In terms of the parameter At$o, which 
is the time interval during which the integrated counts of a burst go from 
5% to 95% of the total integrated counts, the bursts seem to separate cleanly 
into two distinct groups, with the transition occurring around At 90 ~ 2 s. 
This confirms similar indications from earlier experiments (Cline & Desai 
1974, Norris et al. 1984, Dezalay et al. 1992, Hurley 1992), and is the first 
compelling evidence for distinct sub-classes of gamma-ray bursts. (We refer 
here to classical gamma-ray bursts, which represent the bulk of the bursts, 
and do not discuss the class of soft repeaters.) 

An obvious and natural question that arises from this result is: what is 
the relation between the two kinds of bursts? One possibility is that they 
represent two distinct types of sources. For instance, it is even conceivable 
that one class of sources is cosmological and that the other is in the Galactic 
halo. Alternatively, both types of bursts may come from a common source, 
and the differences may arise merely from variations in the initial conditions 
of the sources prior to the bursts, or changes in the environment or the 
viewing angle. We attempt to shed some light on this question by carrying 
out a statistical comparison of the properties of the two kinds of bursts. We 
describe the detector channels on BATSE and some selection effects in §2, 
present the main data analysis in §3, and discuss the implications of the 
results in §4. 
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2 Detector Channels and Selection Effects 



The burst catalog available in the public domain contains a list of all gamma- 
ray bursts that triggered the BATSE detectors between April 1991 and March 
1992 (Fishman et al. 1993). The catalog provides the angular positions for 
260 bursts, the ratio C max /C m i n for 241 bursts, where C max is the maximum 
count rate and C m i n is the detection threshold, the ffuences and peak count 
rates for 260 bursts, and durations for 220 bursts. 

The trigger mechanism and various selection effects of BATSE have been 
explained in detail in Fishman (1992). We will repeat the essentials here. 
The BATSE on-board software tests for bursts by comparing count rates on 
eight large-area detectors to the threshold levels corresponding to three sepa- 
rate time intervals: 64 ms, 256 ms, and 1024 ms. A burst trigger occurs if the 
count rate is above the threshold in two or more detectors simultaneously. 
The thresholds are set by command to a specified number of standard devi- 
ations above the background (nominally 5.5 a), and the average background 
rate is computed every 17 s. Due to an unknown technical reason, many 
bursts have "undetermined" C max / C min in the 256 ms channel. We therefore 
ignore this channel in what follows. 

For a variety of reasons, BATSE has a variable background as a func- 
tion of time, so that the detection threshold C min does not remain constant. 
Moreover, there are periods corresponding to "overwrites" when the sensitiv- 
ity of the detectors is greatly reduced. These effects make an analysis of the 
Cmax/Cmin data somewhat difficult. In order to have a more uniform sample 
to work with, we select a constant threshold C cut , and prune the data so as 
to include only those bursts which satisfy both of the following criteria: 

i r < r 

±. w mln _^ '-'cut; 

2- C max > C cu t. 



The resulting database corresponds to those bursts which would have been 
found by a detector that (i) had a constant threshold of C cut , and (ii) was 
turned on at those times when the real detectors on BATSE had C min < C cut , 
and was turned off whenever the BATSE sensitivity was poorer than C cut . By 
cutting the data in this fashion, we are guaranteed to have a sample of bursts 
with a constant detection threshold and a uniform selection bias. Of course, 
in the process we lose a few bursts, which causes loss of statistical accuracy. 
To minimize this, we select C cut such that the number of bursts retained in 
the database is maximized. This leads to the choice C cut = 71 counts for the 
64 ms channel and C cut = 286 counts for the 1024 ms channel. Fortunately, 
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only a few bursts are eliminated for these choices of C cut ; moreover, the 
excluded bursts are mostly those that are labeled "overwrites", and would 
have been eliminated in any case. 

Following Kouveliotou et al. (1993) we define "short bursts" as having 
At 90 < 2 s, and "long bursts" as having At 90 > 2 s. Before going into the 
main analysis, which we discuss in the next section, we explain first some 
selection effects associated with the different sensitivities of the 64 ms and 
1024 ms channels to the two kinds of bursts. Fig. 1 shows (C max /C cu t)io24 
in the 1024 ms channel vs. (C max /C cu t)64 in the 64 ms channel for the short 
and long bursts. It is apparent that the 1024 ms channel is more sensitive to 
long bursts, while the 64 ms channel is more sensitive to short bursts. Both 
of these effects are quite natural, as we now show. 

The noise in the background counts increases, generally, as the square root 
of the integration time. Therefore the noise is expected to be 4 times greater 
in the 1024 ms channel than in the 64ms channel. Now, if a burst has a 
broad luminosity maximum extending over a time interval greater than 1024 
ms then there will be 16 times more signal counts in the 1024 ms channel than 
in the 64 ms channel, and the signal-to-noise ratio will be 4 times greater. 
Long bursts are likely to display this behavior. On the other hand, if a burst 
is extremely narrow, with a duration less than 64 ms, then the number of 
signal counts will be the same in both channels. The most extreme short 
bursts will correspond to this limit. Based on this argument, we see that the 
ratio (C max / C cut ) 1024 to (C max /C cut )64 for the two channels must satisfy 



with short and long bursts tending towards the lower and upper limit respec- 
tively. 

Comparing the ratio in eq (1) for the long burst bursts, we find 



where the error estimate reflects the width of the distribution. Since R\ 
is not very different from the maximum value of 4, we conclude that, to 
a first approximation, the long bursts tend to have broad nearly constant 
luminosity profiles. For a more detailed description of the intensity statistics, 
we note that the maximum luminosity in the 64 ms channel is greater by a 
factor ~ 4/2.5 = 1.6 than the luminosity in the 1024 channel. This may be 
interpreted as evidence for a "fractal" behavior in the burst luminosity, 




<4, 



(1) 





L(At) oc At 
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where L(At) represents the maximum luminosity of a burst as measured with 
a time constant of At. 

In the case of the short bursts we find 

/ (Cmax/CcutW \ =0 . 76±0 . 48j (4) 

\ { L max /C cut ) 64 / short 

which shows that the 64 ms channel is more sensitive than the 1024 ms chan- 
nel. This is almost entirely because the longer channel dilutes the signal. As 
clear evidence of this effect we note that there is a strong correlation between 
the widths At 90 of the short bursts and the quantity (C max /C cut )io24/(Cmax/Ccut)64 
(the correlation coefficient is 0.49). 

3 Data Analysis 

In this section we carry out a comparison among various samples of bursts. 
We work primarily with three samples: 

• sample S6 4 consisting of 40 short bursts which were detected in the 64 
ms channel, 

• sample Z 1024 consisting of 113 long bursts detected in the 1024 ms chan- 
nel, and 

• sample /io24,64 consisting of 71 long bursts detected in both the 1024 
ms and 64 ms channels. 

In comparing different samples, the key observational data we use are 
the distributions of V/V max of the samples, where V/V max = (C max /C cut )~ 3 / 2 
(Schmidt et al. 1988). Our analysis is based on a simple hypothesis, namely 
that all three populations of bursts have the same underlying spatial distri- 
bution, and that within each population the sources are standard candles. 
On this hypothesis, any apparent differences between the samples are just 
because they are viewed to different distances (or depths) on account of dif- 
ferences in the source luminosity and/or detector sensitivity. Our main result 
is that all the available data are consistent with this idea. 

In our analysis, we make use of the fact that the brighter bursts agree with 
a homogeneous Euclidean distribution, which is characterized by cumulative 
counts varying linearly as V/V max and a mean (V/V max ) equal to 0.5. In 
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contrast, the fainter bursts deviate significantly from such a distribution. 
Based on this observation, we make the following reasonable assumptions: 

(1) We assume that (V/V max ) decreases monotonically with increasing depth 
of a sample. Therefore, if we compare two samples of bursts (with the same 
spatial distribution by hypothesis), we can say that the sample with the 
smaller value of (V/V max ) corresponds to a greater depth or distance. Con- 
versely, if two samples have the same value of (V/V max ) we say that they 
correspond to the same distance. As a further check, when (V/V max ) of two 
samples agree, we compare their full V/V max distributions by means of the 
Kolmogorov-Smirnov test (K-S test), and thus investigate whether or not the 
two populations do indeed have a common spatial distribution. 

(2) If two samples have different V/V max distributions and different mean 
(V/Vmax), but we suspect that they have intrinsically the same spatial dis- 
tribution, then we can bring the two samples into agreement by increasing 
C cut for the deeper population. In effect, we artificially reduce the sensitivity 
of the detector corresponding to the deeper sample so that its sensitivity 
becomes equal to that of the shallower sample. We can think of this opera- 
tion equivalently as reducing the luminosity of all bursts in the deep sample 
by a constant factor. Of course the procedure is meaningful only if it leads 
to agreement in the mean {V/V mauX ) and also in the shapes of the V/V max 
distributions, as discussed in (1) above. 

A point that we would emphasize is that our method of analysis is virtually 
model-free and applies regardless of whether bursts are Galactic or cosmo- 
logical. 

To illustrate the method we consider the two samples of long bursts, Z1024 
and /io24,64- The discussion in §2 showed that the 1024 ms channel is more 
sensitive to long bursts than the 64 ms channel by a factor of 2.5 (eq 2). We 
therefore expect /1024 to correspond to a deeper sample than /io24,64- This 
is confirmed by the mean (V/V max ) values, which are 0.29 ± 0.027 for Z1024 
and 0.38 ±0.034 for /io24,64- Both groups deviate significantly from a uniform 
distribution in Euclidean space, but the deviation is much larger for Zio24- 
We can now artificially bring the two populations to the same distance by 
increasing C cut for the /1024 sample by a factor of 2.5. On doing this we find 
that the value of {V/V max ) for Z 102 4 becomes 0.37±0.035, which is nearly equal 
to the (V/V max ) of the /io24,64 sample, exactly as expected. Furthermore, the 
two cumulative V/V max distributions agree very well with each other after 
this distance correction has been done. The K-S probability (for a worse fit 
than the one obtained) is 86%, which is excellent. These calculations show 
that the /1024 and /io24,64 samples do have the same spatial distribution, with 
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the former being a deeper sample than the latter by a factor of V2.5 in 

luminosity distance. This is no surprise since the two samples have a large 
number of bursts in common, and moreover we know that there is a good 
correlation between the signals in the 1024 ms and 64 ms channels for long 
bursts (cf Fig. 1). The test however demonstrates the validity of the method. 

We next proceed to the more interesting test of comparing the short 
and long bursts. The average (V/V max ) of the s 64 sample is 0.31 ± 0.042, 
as compared to 0.29 ± 0.027 for /1024, which indicates that the two samples 
correspond to nearly the same distance. We now vary C cut of the /1024 sample 
and sq4 sample individually so as to find the range of values over which the 
(VyKiax) values of the two samples agree to within ±la (see Fig. 2). We 
find that 



where the value R3 = 1.4 corresponds to the case when the two (V/V ma , x ) 
values are exactly equal.. As defined here, y/R^ represents the ratio of the 
limiting distances of the /1024 and s 64 samples. We should mention in passing 
that there is a second region of good fit around R3 ~ 6, but the number of 
bursts that survive the cut for this comparison is so small that we do not find 
the solution convincing. Even the primary solution given in eq (5) suffers to 
some extent from the limited number of bursts in the two samples. 

We now test whether or not the S64 and /1024 populations are really consis- 
tent with the same spatial distribution. For this we do a K-S test to compare 
the two distributions of V/V max (see Fig. 2). We see that the test indicates 
fairly convincingly that the two populations do have the same spatial dis- 
tribution. For instance, the K-S probability is 42% if we keep the two C cut 
values unchanged (R% = 1) and 50% if we multiply C cut for the 1024 ms 
channel by the optimum factor of 1.4 to obtain equality of (V/V max ). To 
illustrate the quality of the agreement, we show in Fig. 3 the V/V max distri- 
butions corresponding to the case when R 3 = 1. Note how much better the 
two distributions agree with each other than with the diagonal line which 
represents the homogeneous Euclidean model. The probability that either of 
the observed samples is drawn from the Euclidean distribution, is vanishingly 
small, < 10~ 4 according to the K-S test. 

Finally, we compare the s 64 and /io24,64 samples. In analogy with eq (5) 
we define a corresponding ratio i? 4 for this comparison, 
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For consistency with the two previous comparisons, we expect 



R A — — - — 



(7) 



To check this, we change C cut for the sample from 71 counts to C' cut = 
1.8 x 71 counts and compare the modified s 64 sample with Zio24,64- The mean 
(V/V max ) values of the two samples are 0.41 ± 0.047 for s 64 and 0.38 ± 0.034 
for /io24,64, showing excellent agreement. Further, the K-S test gives a high 
probability of 46%, again in good agreement. We thus conclude that the 
S64 and /io24,64 samples are consistent with a common spatial distribution, 
and that the former is deeper than the latter by about \^L8 in luminosity 
distance. (As an aside we mention that, corresponding to the second solution 
R 3 ~ 6 mentioned earlier, there is an indication of a solution at ~ 0.5, 
but the reduction in counts in these comparisons is somewhat severe and we 
are inclined not to take these solutions seriously.) 

A very interesting feature of the comparison discussed in the previous 
paragraph is that R^ directly represents the ratio of the luminosities of the 
long and the short bursts in the same detector, viz. the 64 ms channel. Since 
the 64 ms channel corresponds to the shortest time interval used in the trig- 
gers, it provides the closest approximation to the instantaneous luminosity 
of a burst. We thus conclude that the maximum instantaneous luminosities 
of the short and long bursts are nearly equal, to within a factor ~ 2. In 
fact, the difference in the luminosities may be even less than the value indi- 
cated in (7). This is because, in the cosmological scenario, the K-correction 
(cf. Piran 1992; Mao & Paczyhski 1992) depends on the spectral index of 
the gamma-ray bursts. Although the spectral indices of bursts have large 
variations, there is some indication that the short bursts are systematically 
harder than the long ones (Dezalay et al. 1992; Kouveliotou et al. 1993). 
The effects of this will be to bring the luminosity ratio even closer to unity. 
Utilizing the fact that we have C min = 71 counts in the 64ms channel, we 
can roughly estimate the peak luminosity of bursts by assuming a power law 
spectrum (see eq 10 below): 



where -D max is the maximum distance to which the bursts can be detected in 
the 64 ms channel. 

An important consequence of the comparisons carried out above is that 
we obtain the relative depths of the short and long burst samples. We can 
therefore estimate the relative number densities in space of the two kinds of 
bursts. We find 



L peak ~ 2 x 10 43 (A„ax/Mpc) 2 erg s~\ 



(8) 




(9) 



8 



We should however mention one caveat, namely that BATSE may have 
missed some very short bursts with durations smaller than 64 ms because of 
dilution. If this is a substantial effect, then the above ratio may be higher, 
possibly closer to unity. 

We now estimate the ratio of the total energy outputs in the short and 
long bursts. We first define an effective duration At cff as the duration of a 
burst would have if it had a constant count rate C max and the same fluence, 
i.e., 

"~ = Wc^ ( ^1JW (10) 

where S is the fluence in units of erg cm -2 in the energy range 50-300 keV 
(which coincides with the trigger energy range), (E) is the mean energy in 
the energy range of 50-300 keV, C max is the maximum count rate in units of 
photons cm~ 2 s _1 , and we assume a power law photon number distribution 
n(E) dE oc E~ a dE. Adopting a = 2 (Schaefer et al. 1992), we find that 
(Atcff) ~ 0.4 s, 12.5 s for the short and long bursts respectively. Combining 
this with the luminosity ratio i?4, we find the total energy ratio is El/Es ~ 
20. It is quite remarkable that the short and long bursts differ by such large 
factors in their durations and total energy outputs, but yet are so similar in 
their maximum luminosities. Finally, for completeness, we mention that the 
short and long bursts are both individually consistent with perfect isotropy. 



4 Discussion 



This paper has been motivated by the recent discovery of Kouveliotou et 
al. (1993) that gamma-ray bursts consist of two distinct subclasses, namely 
short and long bursts. The main aim of our investigation is to test the 
simple hypothesis that both populations of bursts have the same underlying 
space distribution and that within each population the sources are standard 
candles. Our conclusion is that all the available data are consistent with this 
hypothesis. Even though the distributions of C max /C cut for the short and 
long bursts sometimes appear to be different in certain detector channels, we 
are always able to bring the two populations into agreement by modifying 
the detector sensitivity in one or the other sample so as to reduce the two 
samples to the same depth or distance. When we do this, not only do the 
mean (V/V max ) values agree, but also the full distributions of V/V max agree 
when compared by means of the K-S test. Of course, these calculations 
do not prove our basic hypothesis, particularly since we are hampered by 
the small number of bursts in the samples but they make the idea quite 
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plausible. We therefore feel that it is unlikely that one type of bursts (say 
long) is cosmological while the other arises in the halo, or that one is in the 
halo and the other in the disk (Smith & Lamb 1993). It would be too much of 
an accident for the two populations to have the same V/V max distributions. 
Instead, we favor models where the two kinds of bursts arise in the same 
source. In this case, the large difference in durations between the short and 
long bursts may be caused by variations in the initial conditions or in the 
environment of the source or due to differences in the viewing angle. 

Once we accept that the two kinds of bursts arise from a common source 
population, we are able to estimate the relative luminosities and number den- 
sities of the two populations. Surprisingly, we find that both the short and 
long bursts have the same peak luminosity, L pea k ~ 2 x 10 43 (-D max /Mpc) 2 erg s" 
to within a factor of two. The equality of the two peak luminosities is quite 
remarkable when we consider that the durations of the short and long bursts 
differ by ~ 50 and their total energy outputs differ by ~ 20. Incidentally, our 
estimate of L pcak corresponds to the 64 ms channel. On smaller timescales, 
the peak luminosity will be larger, though not by a significant factor (see eq 
3). Further, we estimate the number density of the short bursts to be ~ 0.4 
times that of the long bursts. If the difference between short and long bursts 
is due to viewing angle, this ratio gives an estimate of the relative solid angles 
associated with the two kinds of bursts. 

An important point that we would stress is that these results are obtained 
in an essentially model-independent way. For instance, we do not need to 
make any assumption on whether the sources are at cosmological distances 
or in the Galaxy. 

The constancy of luminosity in two classes of bursts which differ so much 
in their other properties might provide a clue to the physical origin of the 
bursts. Unfortunately, it is virtually impossible to explain the result using the 
most widely used limiting luminosity in astrophysics, namely the Eddington 
limit. If we assume that the sources are not dynamically expanding (but see 
Piran and Shemi, 1993) and take a fixed opacity (e.g. electron scattering), 
the Eddington limit is proportional to the mass of the source, which in turn is 
limited by the variability timescale, St, i.e., L^d < 1-3 x 10 40 (<5t/ms) erg s -1 . 
Taking St ~ 10 ms, the characteristic burst luminosity L pea k that we have 
obtained is marginally consistent with the Eddington luminosity limit for 
sources of mass ~ 1O 3 M located at distances ~ lOOkpc in the Galactic 
halo. Smaller masses are ruled out by the observed isotropy of the bursts, 
while larger masses are ruled out by the variability argument. In fact, Bhat 
et al. (1993) claim to see variability down to 200 /xs, which rules out even 
the 10 3 M Q scenario. 



10 



Another robust limiting luminosity is L max = c 5 /G = 4 x 10 59 erg s -1 , 
which is an absolute upper limit for any source, corresponding to the emission 
of the entire rest mass within a gravitational light crossing time. If we identify 
Lp ea k with this limit, then we obtain the luminosity distance to the sources to 
be ~ 10 8 Mpc, corresponding to a redshift of ~ 10 4 , in an Q = 1 Friedmann 
universe. This is far too large for any known model. 

The fact that it is not easy to come up with a physical explanation for the 
existence of a characteristic peak luminosity for gamma-ray bursts implies 
that, if the effect is real, it may provide an important and vital clue for 
understanding the origin of the bursts. Unfortunately, as we have tried to 
stress, the results have only modest statistical significance at this point. It 
would be very interesting to repeat the analysis with the complete database 
of bursts detected by BATSE. 

We acknowledge helpful discussions with Bohdan Paczyhski. This work 
was supported in part by NASA grant NAG 5-1904. 
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6 Figure Captions 



Fig. 1: C ma x/C m i n for the 1024 ms channel vs. that for the 64 ms channel. 
All bursts with C max / C m i n available in both the 64 ms and 1024 ms channels 
are plotted. Open and filled circles correspond to bursts with durations 
shorter and longer than 2s respectively. Note that C max / C m i n in the 1024 ms 
channel is larger than C max /C m i n in the 64 ms channel for almost all the long 
bursts, while the short bursts tend to have larger C max /C min in the 64 ms 
channel. 

Fig. 2: The thick line shows the probability values Pk-s obtained with the 
Kolmogorov-Smirnov (K-S) test when the /1024 and S64 samples are compared 
versus the quantity logi?3 (see eq 5). Also shown as a thin line is A/ a vs. 
logi?3, where A is defined as the absolute difference of (V/V ma , x ) between 
the two samples, and a is the expected standard deviation. The dotted line 
corresponds to a la deviation. 

Fig. 3: VyVmax for each burst is shown versus the burst's intensity rank 
normalized by the total number. The thick and thin lines correspond to 
the 113 long bursts in the 1024 ms channel (sample /1024) and the 40 short 
bursts in the 64 ms channel (sample s 64 ). The K-S probability that these 
two samples are drawn from the same distribution is 42%. The probabilities 
that these two samples are drawn from a uniform distribution in a Euclidean 
space (indicated as a dashed diagonal line) is < 10~ 4 . 
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This figure "figl-l.png" is available in "png" format from: 
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This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/9304016vl 



This figure "fig3-l.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/9304016vl 



